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Abstract 

This paper deals with the concepts of 7r-flatness and 7r-flat output for hnear 
time-varying delay systems. These notions, introduced and developed by sev- 
eral authors during the last decade, may be, roughly speaking, defined as follows: 
a TT-flat system is a system for which all its variables may be expressed as func- 
tions of a particular output y, a finite number of its successive time derivatives, 
time delays, and predictions, the latter resulting from the advance operator 
TT"^ , TT being a polynomial of 5, the delay operator. Thanks to standard poly- 
nomial algebraic tools, and in particular the Smith- Jacobson decomposition of 
polynomial matrices, we obtain a simple and easily computable characterization 
of TT-flatness in terms of hyper-regularity of the system matrices and deduce a 
constructive algorithm for the computation of 7r-flat outputs. Some examples 
are provided to illustrate the proposed methodology. 
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1. Introduction 



Differential flatness, roughly speaking, means that all the variables of an 
under-determined system of differential equations can be expressed as functions 
of a particular output, called flat output, and a finite number of its successive 
time derivatives ([l^, H, [^, see also 33, 16, 17 1 and the references therein). 

For time-delay systems and more general classes of infinite-dimensional sys- 
tems, extensions of this concept have been proposed and thoroughly discussed in 
lol 24, 31 . In a linear context, relations with the notion of system param- 
28|, |27[ and, in the behavioral approach of ,26], with latent variables 



22 



eterization 

of observable image representations ^35|], have been established. Other theoretic 
approaches have been proposed e.g. in 3^, 3|- Interesting control applications 
of linear time-delay systems may be found in 22|, |2J, [31 1 . 

Characterizing differential flatness and flat outputs has been an active topic 
since the beginning of this theory. The interested reader may find a historical 
perspective of this question in |16l . Il7 |. Constructive algorithms, relying on 
standard computer algebra environments, may be found e.g. in [1] for nonlinear 
finite-dimensional systems, or [2| for linear systems over Ore algebras. 

The results and algorithm proposed in this paper for the characterization 
and computation of 7r-flat outputs for linear time-delay systems are strongly 
related to the algebraic framework developed in [22 , 25 , 31 1 . More precisely, we 
study linear time-delay differential control systems, i.e. linear systems of the 
form Ax = Bu, with x e M" the pseudo-state, and u e M™ the control, for 
given integers m < n, where the entries of the matrices A and B belong to the 
ring A[S, ^] of multivariate polynomials of S, the delay operator, and ^, the 
time derivative operator, over the ground field K of meromorphic functions of 
the variable t. We say that the system Ax = Bu is 7r-flat if, and only if, the 
module generated by the components of x and u over MIS, ^] and satisfying the 
relations Ax = Bu, localized at the powers of a polynomial n G M.[S], is free, 
and a vr-flat output is a basis of this free module (see [22|). 

To characterize and compute 7r-flat outputs, we propose a methodology based 
on standard polynomial algebra, generalizing the one used in [18] for ordinary 
linear differential systems, by extending the original ring M.[S, ^] to the principal 
ideal ring of polynomials of ^ over the fraction field M.{S), namely the 

field generated by fractions of polynomials of S with coefficients in M., and finally 
localize the results of our computations at the powers of a suitable polynomial 
TT of M[S]. This approach allows us to use the well-known Smith- Jacobson (or 
diagonal) decomposition ([^, [l3| ) of matrices with entries in the larger ring 
.^((5)[^] as the main tool to obtain the searched 7r-flat outputs. 

Following , in order to work with a smaller set of equations and variables, 
we eliminate the input variables, leading to an implicit system representation, as 
opposed to previous approaches (see e.g. 22, 24 . Isil. Isl. . Let us also insist on 
the fact that the time-varying dependence of the systems under consideration 
is in the class of meromorphic functions, whereas in jH, HI, this dependence 
is polynomial with respect to time in order to apply effective Grobner bases 
techniques. 
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The main contributions of this paper are (1) the characterization of 7r-flatness 
in terms of the hyper- regularity of the system matrices, (2) yielding an elemen- 
tary algorithm to compute 7r-flat outputs, based on the Smith- Jacobson decom- 
position of the former matrices. In addition, the evaluation of our 7r-flatness 
criterion only relies on computations over the larger ring 

The paper is organized as follows. The 7r-flat output computation problem 
is described in section [2l as well as the algebraic framework. Then, the main 
result of the paper is presented in section [H Finally, the proposed methodology 
is illustrated by some examples in section |4j and its generalization to multiple 
delays is outlined on an example of vibrating string, first solved in (23l |. 



2. Problem Statement 

We consider a linear system governed by the set of time-delay differential 
equations: 

where x G R" is the pseudo-state, u G M™ the input vector, A (resp. B) a 
n X n (resp. n x m) matrix, whose coefficients are multivariate polynomials of 
5 and with ^ the differentiation operator with respect to time and S the 
time-delay operator defined by: 

6: f{t)^Sf{t)^f{t-T), VteR (2) 

where r £ E+ is the delay. 

In order to precise the nature of the coefficients aij{S, ^), i,j = 1, . . . ,n, 
and bij{S, i = 1, . . . ,ri, j = 1, . . . ,rn, of the matrices A{S, ^) and B{S, Jj) 
respectively, some algebraic recalls are needed. 

2.1. Algebraic Framework 

Since we deal with smooth functions of time, a natural field is the differential 
field of meromorphic functions on the real line M. We call this field the ground 
field and we denote it by M.. The previously introduced operators S and ^ 
satisfy the following rules: 

|(a(t)-)=a(t)|+a(i)-, Siaity)^ait-r)5, ^ S-^^ 

for every time function a belonging to .S. The set of multivariate polynomials 
of these operators, namely polynomials of the form 

k,l finite 

is a skew commutative ring [2l|, [s^], denoted by ^[S, The coefficients ajj 
(resp. bij) of the matrix A (resp. B) of system ^ are supposed to belong 
to R[S, thus making system ([ij a linear time- varying time-delay differential 
system, whose coefficients are meromorphic functions with respect to time. 
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2.1.1. System Module, Freeness 

To system ^ is associated the so-called system module, noted A. More 
precisely, following [7, 22], let us consider a non zero, but otherwise arbitrary, 
pair {£,,v) — (^i, . . . , z/i, . . . , fm) and the free modul43j denoted by 
generated by all possible linear combinations of ^ and v with coefficients in 
^[(5, Next, we set 9 = — Bv and construct the submodule [9] of [$,,y\ 
generated by the components of the vector 9. The system module A is, by 
definition, the quotient module A = 

In (22j . in the context of commutative polynomial rings, the notion of projec- 
tive (resp. torsion-free) controllability of a time-invariant system, i.e. a system 
of the form ([T]) with ground field .ft = M, is defined as the projective (resp. 
torsion) freeness of A, and shown to generalize the well-known Kalman control- 
lability criterion to linear time-invariant differential delay systems. Moreover, 
as a consequence of a theorem of Quillen and Suslin, solving a conjecture of 
Serre (see e.g. Il3|), A is free if and only if it is projective free. If F is a 
finite-dimensional presentation matrix of A, the latter module A is projective 
free if F is right-invertible, i.e. there exists a matrix T over ^[5, ^] such that 
FT = I. 

This approach has been generalized to modules over the Weyl algebras by 
Quadrat and Robertz (29I. based on a theorem of Stafford js^l , (see algorithmic 
versions of this result in |l2l[l9j). 

In both time-invariant and time- varying cases, systems whose module is free 
are called flat ( 22, 25, 3lj). Nevertheless, only few systems have a free system 



module, thus motivating the weaker notion of 7r-flatness: we say that the system 
is TT-flat, or that its associated module is 7r-free ([l^, |3l|)jjf, and only if, there 
exists a polynomial tt G ^\S\, called liberation polynomial ([22|,|31|), such that the 
module ^[5, 7r~^, ^] Cg)j^j^ _d_] A, i.e. the set of elements of the form X^ie/ 
with / arbitrary subset of N, G R[5, ^] and £ A for all i ^ I, called the 
system module localized at the powers 0/ tt, is free. In other words, vr-flatness 
means that the state and input can be expressed in terms of the 7r-flat output, 
a finite number of its time derivatives and delays, and advances corresponding 
to powers of the inverse operator it~^ . 

In the sequel, we will also use the extension, as announced, of the ground field 
R to ^{5), the fraction field generated by rational functions of 5 with coefficients 
in ^. The system module over this field extension is ^{5)[-^] A. Indeed, 

freeness (in any sense) of the latter module does not imply freeness (in any 
sense) of the original system module A (see e.g. (2^). 



2.1.2. Polynomial Matrices, Smith- Jacob son Decomposition, Hyper-regularity 

The matrices of size p x q whose entries are in R[6, ^] generate a module 
denoted by Mp^q[S, ^]. The matrix M € Mn,n[S, ^] is said to be invertibl^ if 



^For more details on rings and modules, the reader may refer to [^. 

■^Note that the ^[S, ^]-independance of the n columns and rows of M is not sufficient for 
its invertibility. Its inverse, denoted by A'^, has to be polynomial too. 
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there exists a matrix N G A4n.n[S, ^] such that MN = NM — In, where /„ is 
the identity matrix of order n. The subgroup of A4n,n[S, ^] of invertible matrices 
is called the group of unimodular matrices of size n and is noted Un[S, ^jl- 

Let us give an example of a system of the form that will serve as a 
guideline all along this section to illustrate the various concepts. 



Example 1. 



dt 



where x = {xi,X2)'^, u is scalar and k{t) a meromorphic function. In other 
words ^ reads: 



Xl{t) = k{t){x2{t -T)-X2{t- 2t)) 

X2{t) = u{t - t) 



(4) 



The coefficients o,'>^d —k{t){5 — S^) = —k{t)5{l — 5) are elements oj R[5, ^] 
and the corresponding matrices A and B belong to Ai2,2[5, ^] and M2,i[5, ^] 
respectively. 

Note that it may be necessary to extend the polynomial ring as shown by 
the following computation on the previous example: 

Let us express m of ([3]) as a function of xi. It is straightforward to see 

that X2 = {1 — 5)^^5^^ xi and, since u = 5^^X2, we immediately get 

k dt 

<5-2(i I + ±^ ] !Lx^_ Denoting by tt = {1-5)5'^ G J^[<5], the 




polynomial tt ^ \ ^-r^ + T~r \ ~r li'^'^s in R[5, tt ^,4-1 the ring of polynomials 
\ k dt j dt 

of (5, TT^^ and ^ with coefficients in but not in ^[5, J^]. This is why we may 
also introduce matrices over ^[5,71^^, ^] for some given tt in A[S\. The corre- 
sponding module of matrices of size p y. q will be denoted by Mp^q[5,T:~^ , 
More precisely, a matrix M belongs to J^p q[5, t^~^ ,-^] if and only if there exists 
a finite s e N such that tt" • M e Mp^q[5, ^]. 

Remark 1. It may be argued that the previous expression of u in function of 
Xi is not feasible since (1 — S)^^ = X^T^o ifnp^i^s 



+ 00 



which involves an infinite number of delayed terms. However, if we deal with 
motion planning, if xi is chosen constant outside the interval [^0,^1], for some 



^It is also often denoted by GLn{R[&, -^]) 
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to,ti € M, to < ti, then xi and xi are identically equal to zero before to and 
after ti, and the above series has at most [ *^~*° ] + 2 non vanishing terms, 
where the notation \r~\ stands for the least integer upper bound of an arbitrary 
non negative real r. 

Unfortunately, A[S, ^] and A[S, tt~^, ^] are not Principal Ideal Domains (see 
e.g. 2l|, |32|, [isl, y ) , a property which is essential for our purpose (see the Smith- 



Jacobson decomposition in Appendix |B|. However, if we extend the ground 
field M. to the fraction field A{S), R{S)[-^] is a principal ideal ring of polynomials 
of ^. We then construct the modules A4p,q {S) [^] of matrices of size p y. q 
and Up [5) [^] of unimodular matrices of size p x p respectively, over .S(5)[^]. 
Note that R{5)[j^\ strictly contains J?[5, ^] and J?[(5, 7r"\ ^] for every tt G R[5]. 
Therefore, to interpret the results of computations in ^((5)[^], which turn out to 
be quite simple, and to decide if they belong to a suitable R[5,-k~^ , following 
(22.] , we have recourse to the notion of localization introduced in subsection l2.1.1l 
This aspect will be discussed in section |21 

Since R{5)[-j^] is a Principal Ideal Domain, the matrices of A^p,g {5) [^] enjo" 
the essential property of admitting a so-called Smith-Jacobson decompositio 
or diagonal decomposition: 

Theorem 1 (Smith-Jacobson decomposition [l^). Let M G ■Mp,q (<5) [^] be 
an arbitrary polynomial matrix of size p x q. There exist unimodular matrices 
U €Up (S) [f^] and V eUg (S) [^] such that: 



UMV={ / A, \ .„ (5) 

if p> q. 



(Ap|Op,g-j,) if p<q 

L 

In both cases, e Ma.a {6) [^]; cf = p or q, is a diagonal matrix whose 
diagonal elements (di, . . . , d^, 0, . . . , 0), with s < a, are such that di is a nonzero 
■^-polynomial for i ~ 1, . . . ,s, with coefficients in K{S), and is a divisor of dj 
for all 1 < j < i. 

A constructive algorithm to compute this decomposition may be found in 
Section |B] of the Appendix. 

Given an arbitrary matrix M, we call l-SJ(A/) (resp. r-SJ(A/)), the left (resp. 
right) Smith-Jacobson subset of unimodular matrices U ^Up {5) [^] (resp. V G 
Uq {6) [j^]) such that there exists V lEUq {5) [j^] (resp. U eUp (S) [^]) satisfying 
the decomposition ([5]). 



*we adopt here the names of Smith and Jacobson for the diagonal decomposition to remind 
that it is credited to Smith in the commutative context and Jacobson [l3l[3| for general 

principal ideal domains. 
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Example 2. Consider again system The Smith- Jacob son decomposition of 
B is straightforward: 



U=(^ \ J j , V = l, A, ^6, UBV=(^l 

If we want to eliminate u in (0), using the previous Smith- Jacob son decompo- 
sition of B, we first remark that the second line of U , which will be denoted by 
U2 = (1 0), corresponds to the left projection operator on the kernel of B, i.e. 
U2B — 0. It suffices then to left multiply A by U2 to obtain the implicit form 



F{S, j^)x ^ U2AX = ( A -6))ill ]-U2Bu^0 (6) 



We may also compute a Smith-Jacobson decomposition of F : we first right 
multiply F by i ^ j to shift the 0-th order term in ^ to the left, yielding 



(? l) = {-k6il-S) i) 



and then, right multiplying the result by 

1 

leads to (1 0). The Smith-Jacobson decomposition of F is therefore given by 

UfFVf = (1 0). (7) 



with 

Vf 



l\f-6-\l-6)-^i s-\l-6r^ii 

1 



(8) 



and Uf = 1- 

As previously discussed, the Smith-Jacobson decomposition has been com- 
puted over the ring M.{S) [^] but, according to its result may be expressed 
in the ring A[S, TTp^, ^] with np = (1 — 5)5. 

It is also easy to verify that, since the matrix F is a presentation matrix 
of the system module over M.[S,TTp^, ■^], the latter module, according to is 
isomorphic to any free finitely generated module admitting the matrix (1 0) 
as presentation matrix, which implies that the localized system module is free 
(see Section \ 2.1.1\) . Note also that the polynomial t:f admits non zero roots: 
every T-periodic non zero meromorphic function f of the variable t satisfies 
TTpfit) = f(t — r) — f{t — 2t) = 0. Therefore, the system module is not A[5, 
free. 
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We now introduce a remarkable class of matrices of ^Ap^q {S) [^] called hyper- 
regular. 

Definition 1 (Hyper-regularity [iBl)- Given a matrix M e -^p,? (<5) [^]; we 
say that M is M.{S) [^] -hyper-regular, or simply hyper-regular if the context is 
non ambiguous, if and only if in Ap — Ip if p < q (resp. Ag = Iq if p > q). 

Remark 2. It is not difficult to prove that a finitely generated module A over the 
ring A{S)[-^] whose presentation matrix F is hyper-regular cannot have torsion 
elements and is therefore free, since ^((5)[^] is a Principal Ideal Ring. However, 
the system module A, over ^[6, J^], does not need to he free, as shown in the 
previous example. Nevertheless, it will he seen later that there exists a liberation 
polynomial it, deduced from the Smith- Jacohson decomposition of F , such that 
the ^[(5, TT^^, -j^]-system module, associated to the same presentation matrix F, 
is free. 

Example 3. Going back to the decomposition of Example [H a variant of this 
decomposition may he obtained as: 

J ), V^l, A;=/i = 1 (9) 

which proves that B is ^{5) \^-^ -hyper- regular. 

It is also immediate from ^ that F is ^{5) [-^] -hyper-regular, or more 
precisely, ^ [6, iTp^, ^] -hyper-regular with np ~ {1 — S)S. Note, on the contrary, 
that A is not K{5) [^] -hyper-regular, though F — U2A is. It is easily verified 
that 

UaAVa = Aa 

with 

^ i-l'+i kS{l-5)) 

and 

" \ f-^ + J ■ 

\ " k ^ dt) dt / 

Thus, since the second degree -^-polynomial (^f + of the diagonal of A a 

cannot be reduced, A is not hyper-regular. 

2.1.3. Implicit system representation 

One of the applications of the Smith- Jacobson decomposition concerns the 
possibility of expressing the system ([1]) in implicit form by eliminating the input 
It, which may be useful to work with a smaller number of variables. 

For simplicity's sake, we rewrite system ([1]) Ax = Bu. 
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Proposition 1. System |ip is equivalent to 

FX^O, AbN-^U^ {I,r^,Om,n-m)MAx (10) 

with 

In^.m)MA, (11) 

M e l-SJ(B) and N such that 

MBN = ( ) . (12) 

y Un— m,m J 

Moreover, if B is ^{5) [^] -hyper-regular, the explicit form flj) admits the 
implicit representation 

Fx^O (13) 

with F given by ill]) , and with Ab — Im in ilS^) . In this case, u is deduced 
from x by 

U = N{Ira,Om,n^m)MAx. (14) 

Proof. Consider a pair of matrices M and N obtained from the Smith- Jacobson 
decomposition of B, i.e. satisfying (|12p . Thus, left-multiplying both sides of 
system H]) by {On-rn,m, In-rn)M , according to ^U} we get Fx — 0. On the 
other hand, multiplying both sides of ^ by {Im,Om.n-m)M we get 

{I^,Om,n-rn)MAx = AbA^~^U, 



hence the representation ([TI 

Conversely, if x and u are given by (jlOp . we have 



MAx=[ ^'^^^^ ).= ( ^^"f"" ]=MBu 



(On — m.m 1 -^n— m)-^^^ / V ^ 

the last equality being a consequence of (|12p . Thus, since M is unimodular, the 
pair (x, m) satisfies Ax = Bu, which proves the equivalence. 

Finally, if B is M.{S) [^] -hyper-regular, one can replace by and the 
second equation of (fTU|) becomes (dH). Thus, u is a M.{S) [^]-combination of the 
components of x and can be eliminated. Therefore, the remaining part (1131) is 
the desired implicit representation of ([1]) . The proposition is proven. □ 



In the sequel, if B is hyper-regular, we refer to (jl3p as the implicit represen- 
tation of system ([T]) . 

Proposition 2. For system {Ip to be M.{d) [-j^] -torsion free controllable (see 
^1; El/J; it is necessary that B and F, defined by \11\) . are ^(S) [■^'^-hyper- 
regular. 

Moreover, in this case, there exists a polynomial 7f such that the localized 
system module at the powers of 7f is free. 
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Proof. Assume that the system is A{S) [^] torsion free controhable and that B 
is not A{6) [^] -hyper-regular. Then, using the decomposition has at 

least one diagonal element which is a polynomial of degree larger than or equal 
to 1 with respect to ^ and with coefficients in R{S). There indeed exists a non 
zero element v of the R[6, ^]-free module generated by the components of u, 
such that Abv = (v is SL non trivial solution of a differential delay equation). 
It is immediately seen that the pair x — 0,u — Nv is a non zero torsion element 
of the system module, which contradicts the freeness assumption. 

If F is not Si{S) [^]-hyp er-regular, its decomposition is given by UFQ = 
{Ap,0), Ap having at least one diagonal element which is a polynomial of degree 
larger than or equal to 1 with respect to ^ and with coefficients in A{S) which 
shows, using the representation pH)) . that every pair (^i, ^2) such that Ap^i = 0, 

^1 7^ 0, and ^2 arbitrary, satisfies (Ap,0) ( ) = 0. Thus, the pair (x,u) 
with X = q( I and u satisfying AbN^^u — (Ira,Om n-m)MAQ ( 

V « y ' V « 

(see ((TU)) ) is a torsion element of the system module. Consequently, the latter 
module cannot be R{S) [;^]-free. 

To prove the existence of tt , we remark that, according to the Smith- Jacobson 
decomposition algorithm (see Section |B] of the Appendix), if M £ \-5^{B) and 

N £ r-SJ(i3) are such that MBN = ^ ^q" ^ : each row of M and N may con- 
tain the inverse of a polynomial of ^[S\. Taking the LCM, say TrM,N, of these 
polynomials for all rows, we immediately get that ttm.n • M G W„[(5, 'Km,n • 
N e Um[5, ^] with 'n:M,N G ^[S\- The same argument applies for F: con- 
sider U e Wn-rn (i^) [^], Q G Z//„ ((5) [^] as abovc, namely such that UFQ = 
). There exists tTjj q e J?[(5] such that Tr^ g • U G W„_,„[(5, ^], 

TTjj Q ■ Q £ Un[S, J^]. Taking tt as the LCM of ttm.at and t^jj q, it is immediately 
seen from the decompositions of B and F, multiplied by suitable powers of n, 
that the system module over the localized ring M.[S, , ^] is free, and the proof 
is complete. □ 

2.2. Differential n -Flatness 



We first recall the classical definition of a flat system [16|,|33|, in the context of 
systems described by ordinary nonlinear differential equations: a system is said 
to be differentially flat if and only if there exists a set of independent variables, 
referred to as a flat output, such that every system variable (including the input 
variables) is a function of the flat output and a flnite number of its successive 
time derivatives. More precisely, the system 

X = f{x,u) 

with X E R" and u E M™ is differentially flat if and only if there exist a set of 
independent variables (flat output) 

y = h{x,u,u,u,...,u^'''>), yE'MT (15) 
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such that 

^ = oi{y,y,y,...,y^'''>) (16) 
u = /3(y,2/,y,...,y(^+i)) (17) 

and such that the system equations 

^(y, y, j/, . . . , 2/(^+1)) = / [a{y, y,y,..., y^'^),f3{y, y,y,..., y^'+'^)) (18) 

are identicahy satisfied for ah smooth enough function t t-^ y{t), r and s being 
suitable finite m-tuples of integers. 

Let us now discuss further extensions of this definition to hnear delay systems 
and consider first an elementary example to motivate the next definition. 

Example 4. Let us consider the elementary system 

x{t) = u{t — t) 

with A = 2i ^i^d B = 5 in the notations of Clearly, if we set y — x, 

y looks like a flat output though u = 5^^y contains a one-step prediction and 
belongs to R{5)[-^] but not to ^[5, J^]. However, for motion planning, such a 
dependence remains acceptable (see Remark]^, even if for feedback design it 
poses more delicate problems. This notion is called differential S -flatness (see 
e.g. MM])- 

According to [2^, this notion is generalized as follows: 

Definition 2 (Difl[erential 7r-flatness [l^l). The linear delay system (QP is said 
to be differentially n-flat (or n-free) if and only if there exists a polynomial 
TT G ^[S], and a collection y of m {5,TT~^)-differentially independent variable^, 
called n-flat output, of the form 

y^Poid,7r-')x + PiiS,7r-\j^)u, (19) 

with PQ{S,Tr^^) G (.^[(5, vr^"'^])™^" the set of matrices of size m x n, with coeffi- 
cients in A[6, TT^"'"] , Pi{S, TT^^ , ^) e M.m.m[S, TT~^ , and such that 

X = Q{5,n-\j^)y, (20) 
u = R{5,n-\j^)y, (21) 

with Q{d,TT-^,f^) G Mn,m[S,Tr-\f^], and R{S , TT''^ , f^) E Mm,,n[S,n-\ f^]. 



^more precisely, there does not exist a non zero matrix S £ Mm,,m[S,'rv ^, ^] such that 
Sy = 0, or equivalently, Sy = impUes 5 = 
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In other words, definition [H states that the components of a 7r-flat output y 
can be obtained as a ^[(5, 7r~^, ^]-linear combination of the system variables, and 
that the system variables {x,u) are also ^[(5, tt"^, ^]-linear combinations of the 
components of y. Thus x and u can be calculated from y using differentiations, 
delays, and predictions (coming from the inverse of tt). Note that the facts that 
every element of the system module is a ^[5, 7r~^, ^]-linear combination of the 
components of y, and that the components of y are ^[(5, tt"-'^, ^] -independent, 
indeed imply that the components of y form a basis of the system module 
R[5,'K~^, ^] A.] A, which is therefore free, hence the equivalence with the 
definition of subsection 12 .1.1 1 

If system ([1]) is considered in implicit form (|13p after elimination of the 
input M, since this elimination expresses u as a J? [5) [^]-combination of the 
components of x, the expression ([T^. combined with P^]). reads y — Px with 
P G Mm.n [5) [^]- The previous definition is thus adapted as follows: 

Definition 3. The implicit linear delay system hlS^) is said to he differentially 
TT-flat if and only if there exists a polynomial tt e A[S], and a collection y of m 
{6,TT^^)-differentially independent variables, called n-flat output, of the form 

2/ = F(<5,7r-\|)x, (22) 

with P{6,n^^, ^) e A4m,n[^,T^^^ , and such that 

x = Q[5,Ti-\^)y, (23) 

WlthQ{5,TT-\j^)eMnA5,^'\ji]- 

The matrices P, Q and R of p^ - (PT|) in the explicit case, and P and Q of 
(|22p - (j23p in the implicit case, are called defining operators of the vr-flat output 

y- 

Remark 3. In f2^] and later (see e.g. fslij), the above notion is often called 
■n-freeness and introduced via the notion of system module. The wording tt- 
flatness appears, to the authors knowledge, for the first time in [24]. It has 
also been related to system parameterization in Q/. We have preferred here 
the name n-flatness, in reference to differential flatness, and to directly present 
it via the notion of flat output, rather than basis of the system module. Note 
that in formula filS^) . Po is a 0th degree polynomial of ^ to mimic the general 
definition in U5]} that does not include time derivatives of x, with il9\} ~ ^21]) 
restricted to linear expressions. 

Example 5. Let us go back again to example 1 and let us prove that y ~ xi is 
a TT-flat output with tt = (1 — S)S'^. From we have 

X2 = d-\l - 6)-'jXi = d-\l - 6)-'jy 



12 



In other words, following the notations of (P3'-([1IP, Po = (IjO), Pi = awrf 



u = n ^ 



k d l(i'^\A„/,- iC?N 



which proves that y is a 'K-flat output. 



3. Main Result 

In this section, we propose a simple and effective algorithm for the com- 
putation of TT-fiat outputs of linear time-delay systems based on the following 
necessary and sufhcient condition for the existence of defining operators of a 
TT-flat output. Moreover, explicit expressions of tt and of these operators are 
obtained. 

Theorem 2. A necessary and sufficient condition for system (Qp to he n-fiat is 
that the matrices B and F are A{S) [^] -hyper-regular. 

We construct the operators P , Q and R and the polynomial tt as follows. 

0. According to Propositions]^ and\^ construct the decomposition of B ilS^) . 
define F by Ul\) and compute n; 

1. g = Q ( ] , with Q e r-SJ(P). Note that n ■ Q e Mn,m[S, ^] by 
construction; 

2. i? = N{Im,^rn,n-rn)MAQ, with N £ r-SJ(B). There exists ttr G A[S] 
such that TTfl • i? s A4„i^rn[S, 

3. P = W{I,n,Om.n-m)P, With P G l-SJ(Q) and W G r-SJ(Q). There exists 

Tip e R[5] such that np ■ P £ Mrn,n[^, 

Let TT be given by tt = LCM(7f, ttp, tt^), the least common multiple of tt, ttp and 

Thus, P, Q and R are defining operators with y — Px, x — Qy, u — Ry, 
and y is a TT-flat output. 

Proof. If i? or P are not ^{S) [^] -hyper-regular, according to Proposition^ the 
R{S)[-^] system module cannot be torsion free. Therefore, system ^ cannot 
be TT-flat, for any tt G ^[S]. Taking the contrapositive of this statement, the 
hyper-regularity of B and P is proven to be necessary. 
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We now prove that the A{S) [^] -hyper-regularity of B and F is sufficient to 
construct the defining matrices P, Q and R as well as a liberation polynomial 

TT. 

We first use the implicit form (|13p of Proposition [1] and more precisely 
(fTT |) -(fT4 |) to obtain a decomposition of B and F. Since F is hyper-regular by 
assumption, there exist U € Un-m (S) G (^) [^] such that UFQ = 

{In-m On-rn,m)- Consequently Q = Q ^ On-m^m ^ ^ ^.^j^ g ^ |._SJ(^)^ jg such 

that FQ — On-m.m- Thus setting x — Qy we have that FQy — for all y. 

The existence of tt e K[5] is proved following the same argument as in the 
proof of Proposition [2j for an arbitrary a x /3 hyper-regular matrix M, accord- 
ing to the Smith- Jacobson decomposition algorithm (see Section [B] of the Ap- 
pendix), if Um e l-SJ(Af) and Vm e r-SJ(A/) are such that UmMVm = (/q,0) 

if a < /3 (resp. UmMVm = ^ ^ ^ if a > /3), each row of Um and Vm may 

contain the inverse of a polynomial of ^[6]. Taking the LCM, say ttm, of these 
polynomials for all rows, we immediately get that ttm ■ Um & ■M.a,i3[S, with 
ttm S ^[S]. Applying this result to the decompositions of B and F, we have 
proven the existence of tt such that n ■ Q £ h{n[S , and if ■ Q £ Mn,m[S, ■^], 
which proves item 1. 

Going back to and p^ . setting R ~ N {I,n 0„i,n-m) MAQ, we obtain 
u — Ry with N G r-SJ(i?). Finally, the proof of the existence of ttj^ G ^[5] such 
that nji ■ R £ M.,n,m.[S, ^] follows the same lines as in item 1, which proves 2. 

Since Q is hyper-regular by construction, its Smith- Jacobson decomposition 
yields the existence of F G W„ (6) [^] and W G Um (S) [ Jj] such that 



PQW = 




thus {Im Om,n-m) PQW = Im and. Setting 

P=W{Im,Qrn,n-rn)P 

it results that Px = W{Im, Om,n^m)Px = W{Im,Om,n-m)PQy = WW-'^y = y, 
and the proof of the third item is complete, noting again that the existence of 
TTp G M.[S] such that np ■ P £ Mm,n[S, ^] follows the same lines as in items 1 
and 2. 

Finally, taking tt = LCM(7f, Trp, ttr), the least common multiple of tt, ttp and 
TTR, it is straightforward to show that n ■ P £ Mm,n[S, n ■ Q £ Mn,m[S, ^] 
and TT-R £ Mm,m[S, Therefore, P, Q and R are defining matrices of a 7r-flat 
output for system ([1]) and thus that system ([T]) is 7r-flat. □ 

Theorem [5] is easily translated into the Algorithm [T] presented below. 

Remark 4. The -K-flatness criterion is given in terms of properties of the ma- 
trices B and F that depend only on the larger ring ^{6)[-^]. If, in addition, 



14 



Algorithm 1: Procedure to compute 7r-flat outputs. 

Input: Two matrices A e Mn,n[S, ^] and B e Mn,m[S, 
Output: a polynomial tt e ^[(5] and defining operators 

R G A^m,m['^, TJ""^, ^] such that y = Px, x — Qy and u — Ry. 
Initialization: Test of hyper-regularity of B by its Smith- Jacobson 
decomposition, which also provides M G l-SJ(i3) and 
N e r-SJ(B), i.e. such that MBN = , On-m,mf ■ If 
B is not hyper-regular, the system is not 7r-flat whatever 

Algorithm: 

1 Set F = {On-m,m, In-m)MA and test if F is hyper-regular by computing 

its Smith- Jacobson decomposition: VFQ = (J„_m, 0„_m,m)- If F is not 
hyper-regular, the system is not 7r-flat whatever tt £ ^[5]. Otherwise, 
compute 7f e j?[(5] such that tt • M e Z^„[(5, tt • TV e W„, [(5, 

2 Compute Q = Q (0„_„,™ , /m)^, i? = A^(/™, 0„,„_„)AfA(5 and 

TTB. e J?[(5] such that ttr ■ R e Mm,,m[S, 

3 Compute a Smith- Jacobson decomposition of Q: 

PQW = {Im , 0„_™,m)^, P = Ty(Jm, 0™,„_m)F, and a polynomial 
vrp e J?[(5] such that ttp • F e Xm,n['5, 

4 Compute the polynomial vr = LCM(7f, ttp, tt/?) G .R[(5]. The system is 
TT-flat. 



the suhmodule of ^[S\ generated by the powers of tt is torsion free (e.g. t: — 6, 
for which the equation Sf = admits the unique solution f = 0), then the sys- 
tem module A over A[S, tt"^, ^] is torsion free, and the original module A (over 
^[6, -^]) is free if, and only if tt — 1. Note that computing a free basis of the 
system module directly over M.[S, -^], would require more elaborated tools such as 
those developed in fM, H [ID- 
Remark 5. Since the computations are made in the larger modules Aip^q {5) [^] 
for suitable p and q, Theorem [H remains valid for systems depending on an 
arbitrary but finite number of delays, say 5i, ... ,5s, by replacing the field R{5), 
on which Mp,q (S) [^] is modeled, by the fraction field M.{6i, . . . , 6s) generated 
by the multivariate polynomials of di, . . . ,ds. An example with two independent 
delays is presented in Example \4-.S\ 

Remark 6. Our results also apply in the particular case of linear time-varying 
systems without delays. If H is the field of meromorphic functions of time, 
is a Principal Ideal Domain and all the computations involved in Theorem\^and 
the associated algorithm, remain in this ring, contrarily to the case with delay. 
Therefore, the last step consisting in finding the so-called liberation polynomial 
TT is needless. Related results may be found in Jl^ . [^ /. 
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4. Examples 

Back to the Introductory Example 

Going back to Example [U let us apply the previous algorithm to the time- 
delay system defined by ([3]), for which a 7r-flat output is already known from 
Example [5] The first step, consisting in the computation of a Smith- Jacobson 
decomposition of the matrix B has already been done, the left and right unimod- 
ular matrices M e l-SJ(B) and N G r-SJ(B) being given by dH) with M = U' 
and N = V ~ I. Then the matrix F = {0 1)MA of an implicit representa- 
tion of is given by ([5]). Its Smith- Jacobson decomposition VFQ = (1 0) 
is given by ([T])-®, with V — Up — I and Q = Vp, and has been seen to be 
hyper-regular in Example [3] We easily check that tt ~ (1 — &)5. 

Going on with step 2, we set 



and we have -kb. — t^S. 

Note that, setting x — Qy and u = Ry, we recover formulae ([25]) and (|24p 
(or equivalently (1^ ). 

According to step 3 of the algorithm, we compute a Smith- Jacobson decom- 



hence P = (1 0) and y — Px — xi. Here, np — 1. Finally, the least common 
multiple of (1 — 6)6, (1 — 6)6'^ and 1 is vr = (1 — 6)6'^. We have thus verified that 
Algorithm [T] comes up with the same conclusion as Example [S] 

4-2. A Multi-input Example 

Let us consider the following academic example of multi-input delay system: 



' xi{t) + ' (t) - 2xf ^ {t - r) + ^ (t) + x*"*' {t - t) - (t) + a;*"' (t) - (t) 

-X4, {t) + {t)=Ui{t)+Ui (t) +U2{t), 

ii (t) + ii (t - r) - ii (t - 2r) + a;f ' (t) + ' (t - r) + {t - 2t) - {t - r) 

+2X2{t) +X2{t-T)- xf\t) - xf\t) + ±3(t) + xf\t - r) - XA[t) ~- Xi(t - t) 

-xf'\t)^u^^\t-T)+U2{t-r), 

~xi{t - 2t) + xi{t~ 3r) + (t - 2r) - a;2(t - r) + ±2(4 - 2r) + 2;^^'(t - r) 
-2;3(t - r) + i3(t - 2r) + X4{t ~ t) = Ui{t - 2r) + U2{t - 2r), 

, Xi(t - r) + X2{t) + X-i{t) = Ui{t) + U2{t). 




and verify that tt • Q e A^2,i['^: Moreover 






(27) 
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Denoting by x the state vector, x = {xi,X2,X3,X4)'^ , by u the input vector, 
u = {ui,U2)^, and by S the delay operator of length r, system (|27)) can be 
rewritten in matrix form Ax = Bu, with A G M4a[S, ^] and B G M.i^2[S, ^] 
defined by 

^ = ( A2 A3 At) (28) 



with 



Ai 



A. 



A, = 



and 



A 

dt 



1.(1-25) + 
^(l + <5-<52) + ^(l 

4-J 



IF" 



dt^ 



dt 
d^ 

IP 



dt 



dF 



d_ 

dt 



dt 



dt^ 



dt^ 
_ dx2 



dt^ 



^ dt ' 



B 



A4 



1 + — 

dt 

dt" 
A 
dt 



\ dt 



dt* 



{1 + 6) 



-2-5 
dt" 



dtp 



1 

dt" 
52 



\ it 1 / 



(29) 



We apply Algorithm [T] to compute a 7r-flat output if it exists. We start with 
the Smith- Jacobson decomposition of B. By left multiplying B by the following 
product of unimodular matrices M — M4M3M2M1 e Ui (6) [^], given by 



/I -1 \ 

10 

10 

\ 1 / 



M2 = 



1 



\ 

1 



dt^" 

-is' 1 
V 1/ 



^ 1 ^ 

1 

10 

V 1 y 



( 1 




V 





1 

-it^ 1 

-52 



\ 



1/ 



and by setting N = I2, we obtain the Smith- Jacobson decomposition 



/ 1 



MBN = 







1 
1 



dt "* ^ 
dt" 



1 

a 
dt 



1 













dt" 









1 




52 














1 


) 




I 





) 
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thus showing that B is hyper-regular. We then compute an imphcit represen- 
tation of (EZl) by = (02,2 I2)MA, i.e. 



F — { Fi F2 F3 F4 



(30) 



with 



F3 



dt 



dt'i 



dt^ 



(31) 



d_ 

dt 



F4 = 



is 



to which corresponds the diiTerence-differential system 

(xi {t) +xi{t~T)- xi it - 2t) + xi {t) + £ 1 (t - r) - x^^' (t - r)) 

+ (2X2(0 +X2{t-T)- i2{t) -X2it-T)- X^^\t)^ 
+X3{t) - (X^it) + Xi{t - t) + Xi{t)) = 0, 



i-xi{t - 2r) +xi{t- 2t)) - (.X2(i - r) 



•^2 



-.T3(f - t) + X4{t - t) = 0. 

According to step 1, we compute a right Smith- Jacobson decomposition of F: 



VFQ 



where V = 1 and 

/ 




10 
10 



(32) 



V -(1 + '^)" 




1 

4-1 

dt^ 

dt) dt 



d-' 
dt^ 

(A 
\dt 



1 



dp 



showing thus that F is hyper- regular and that tt = d{l + 6). 

For the interested reader, Q is obtained as the product Q 
matrices of elementary actions: 



Q1Q2Q3 of 



Q2 







/ 





1 

















1 








1 - 




dt ^ 







•} 








V 


1 





} 








91,1:2 


91,2,2 


91,3,2 




91.4, 2 


\ 









1 




















1 









V 















1 


/ 



18 



with 



91,1,2 
9l,2,2 

91 



,3,2 = (l + <5)-^ ((2 + '^)l-^)-| 



and 



9l,4,2 



Q3 



dt " 



I 1 



V 



(1 + -^) 







-1 _d |_ u 



According to step 2, we get 



Q = Q 



02,2 

/2 



-1 + 
1 





V (1 



dF 



dt J 



dt- 



dt 






1 



1 



d'' I d^ 
dF ' dp 

Ut ^ 



dF 



-<5 



(33) 



i? = iV(/2 02,2)MAQ : 



dF 



dF 



d_ _ d-' _ d; 
„ dt dt'^ dt 

d^ , d=* I 2 d^ 

"dF dF d? 



+ 



dV> 



with TTfl 



and 



Ulit) 



1. Prom X 
( 



Qy and u = Ry, we deduce the expressions 



y2{t) 

yiit) 

yi (t) + ^ it) - 2/2 (i - r) + y (t) + yf^ (t) 

V yi{t)-y?Ht) + y2it)-y2{t-r) + yPit) J 



-y['\t) + Ut)-yi'\t)-yi'\t) 

y['\t) + y['\t) - y^^\t) + yi'\t) + 2yi'\t) + yi'\t) 



(34) 



(35) 

Next, according to step 3, we compute P € U4 (S) [4t] and W G U2 (S) [-i] 



such that PQW = [h 
( 



V 



_d_ 

dt 
dS 
dt 



d^ 
dp 



S - 



dp 



d_ 

dt 



d" 
dp 
_d 
dt 



d^ 
dp 



\ 

* 

1 
1/ 



W = h. 



Again, P is obtained as the product P = P2P1 of elementary actions: 



P2 = 



V -(i-i)i -d+i-^) 



dt 



\ 



1 
1 / 



/ 1 \ 

10 

10 

\ 1 / 
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Then 

P = W{l2 02,2) P = 
which, with y = Px, yields 



10 
10 



2/1 = X2, J/2 = Xi 



(36) 



and TTp = 1. 

Then it is immediately seen that tt = (5(1 + 5) and that the system is 7r-flat. 

Remark 7. It is worth noting that the polynomial tt = 6{1 + S) only appears at 
the intermediate level of the computation ofQ, and not anymore in the defining 
matrices P , Q and R. However, this means that the system module contains 
elements z such that z{t) = ~z(t ~~ t) /or all t, that satisfy ttz = 0, thus 
preventing this module from being free. 

4.3. Vibrating String With an Interior Mass 

As noted in Remark [5l the computation of 7r-flat outputs based on Theorem 
[5] can be extended to linear systems with multiple delays. As an example, 
we consider the system of vibrating string with two controls proposed in |23| . 
which can be modeled as a set of one-dimensional wave equations together with 
a second order linear ordinary differential equation describing the motion of 
the mass. Using Mikusihski operational calculus (see for instance ^1^), this 
infinite-dimensional system can be transformed into the time-delay system 



Mt - 2ti) + iPiit) ^ Ul{t - Tl), 
Mt) + ^2{t - 2t2) ^ U2it ~ T2). 



(37) 



where ryi and 772 are constant parameters. Denoting the state a; = (V'l, (^i, '02, '^2)'^, 
the control input u = (ui,U2), and Si, 82 the delay operators of respective 
lengths Ti and r2, the system p7p may be rewritten in the form Ax = Bu, with 
A e 7^4,4 {5i,S2) and B e M^a {Si, 82) [^] given by 



A = 



( . ' 


1 


-1 


-1 






( ' 





\ 


i + m 




'72 


-m 














1 



















V 





S'2 


1 


) 




I 


S2 


J 



(38) 



The computation of a Smith- Jacobson decomposition of B is here straightfor- 
ward: it suffices to exchange the two last lines of B with the two first lines, and 
h 
02,2 



we get MBN 



with 



M 



/ 








1 















1 















V 





s-^ 





/ 



N = l2. 



(39) 
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We have then ttm.at = 6162- Thus B is hyper-regular and 



F = (02,2 I2)MA = 



A right Smith- Jacobson decomposition of F, namely VFQ = {I2, O2 2), is given 

by 



Q = 



v = 



Si 

001 
\ 1 / 



(40) 



with TT = ^1^2, and where Q is obtained as the product of elementary actions 
Q\ and Q2: 



( 1 








1 


2jji 








V 







(41) 



thus showing that F is hyper-regular. 

According to step 2 of the algorithm, we compute Q 
R = N{l2, 02,2) MAQ: 



02,2 
I2 



and 



Q = 



2i(^ + (^l+%)) A{i + ivi-V2)) 

1 

V 1 J 

SI 1 



(42) 



R 



with 



= ^i-i + ivi + %)) + + m)) 

We indeed have TTij = 1. 



(43) 
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Therefore, setting x — Qy and u ~ Ry, we obtain the expressions 





( 


Mt) 








V Mt) J 


V 



- V2)yi{t) - hit) 
m)vi{t) + hit) ^ 

2/1 W 
2/2 W 



- {til +m)v2{t)) \ 
ivi - m)y2{t)) 

) 



(44) 



and 



= 2^ [2/1 ^ "^l) ^ 2/1 + '^1) + 2/2(^ - Tl) 

Hvi + m){yi{t - n) + y2{t + n)) + (771 

U2{t) =yi{t- T2) + y2{t + T2). 



hit + Ti) 

- mKviit + ri) + y2{t - n))] 

(45) 

Further, according to step 3, we compute P and of a Smith- Jacobson 

I2 



decomposition of Q, namely PQW — 
hhhPi with 



Pi 



( 

1 

V 



( 1 





V 1 






1 / 

1 


/ 



P2 



02 



1 

2>)i 



We find W ^ I2 and P = 



'A. 



thus 



P = 



/ 



1 

Vo 1 



1 

2r;i 
^2^ 



\dt 



Pi 



1 



-(^1 



\ 

I 1 



V 



^2)) 

^'72)) 



2))i \ dt 



1 

+ (??1 

- (?7i - 







^??2)) 
'72)) 







1 

1 








2^(i-(^H 
('71 



1 

"2))i 



and 



P = T4^(/2 02,2) P = 



1 M _ 

2))i Vdt 
1 M 
2r)i Vdi 

3 1 
3 




1 

- ivi 

+ {Vi 


1 





1 / 



1 / 



^2)) 

-V2)) J 



(46) 



(47) 



Finally, setting y — Px, we get yi — 1^2, y2 — 4>2 and ttp = 1. 
Taking the least common multiple of (1, 1, 6182), we get vr = 5i52- It is then 
immediately seen that yi = ip2 and j/2 = 02 is a (5i(52-flat output. 

Remark 8. For multiple delays 61,62,. ■■ , 6n, a flat output for which the poly- 
nomial IT is restricted to a monomial (i.e. of the form ^^^^2^ • • '6^") is called 
6- flat output in Wii l. This is the case here with n = 2 and si = S2 = 1. 

In fil /. a different solution yi = 6i(j)i — w, 2/2 = '/'i + ^1 has been 



Remark 9 

proposed. 
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5. Concluding Remarks 



In this paper, a direct characterization of 7r-flat outputs for hnear time- 
varying, time-delay systems, with coefficients that are nieromorphic functions 
of time is obtained, yielding a constructive algorithm for their computation. The 
proposed approach is based on the Smith- Jacobson decomposition of a polyno- 
mial matrix over the Principal Ideal Domain ^{S)[-^], containing the original 
ring of multivariate polynomials K[S, The fact that the computations are 
done in a larger ring, which is a principal ideal ring, makes them elementary 
and their localization at the powers of a (5-polynomial results from an easy 
calculation of least common multiple of a finite set of polynomials oi S. It is 
remarkable, however, that the 7r-flatness criterion only involves properties of 
the system matrices over the extended ring Several examples are pre- 

sented to illustrate the simplicity of the approach. Translating our algorithm 
in a computer algebra programme, e.g. in Maple or Mathematica, might be 
relatively easy and will be the subject of future works. 
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Appendix 

We recall from section [OI the following theorem: 

Theorem 3 (Smith- Jacobson decomposition [^, [111), -^^i M G Mp.q {5) [^] be 
an arbitrary polynomial matrix of size p x q. There exist unimodular matrices 
U €Up ((5) [|^] and V GUq (S) [^] such that: 

• UMV^{Ap\Op,g_p) tfp<q, 

. UMV= ( ) tfp>q. 

In both cases, Ap € Mp,p {S) [^] and Aq € M.q,q {5) [^] are diagonal matrices 
whose diagonal elements (c?i, . . . , do-, 0, . . . , 0) are such that di is a nonzero 
polynomial for i — 1, . . . , cr, with coefficients in ^{6), and is a divisor of dj for 
all i < j < (7. 
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A. Elementctry Actions and Unimodular Matrices 

The group of unimodular matrices admits a finite set of generators corre- 
sponding to the following elementary right and left actions: 

• right actions consist in permuting two columns, right multiplying a column 
by a non zero function of ^{6) , or adding the j th column right multiplied 
by an arbitrary element of R{S)[-^] to the ith column, for arbitrary i and 

i; 

• left actions consist, symmetrically, in permuting two lines, left multiplying 
a line by a non zero function of ^((5), or adding the jth line left multiplied 
by an arbitrary element of ^{S)[-^] to the ith line, for arbitrary i and j. 

Every elementary action may be represented by an elementary unimodular ma- 
trix of the form Tij(j}) = + li.jP with lij the matrix made of a single 1 at 
the intersection of line i and column j, 1 < i,j < v, and zeros elsewhere, with 
p an arbitrary clement of A{d)[-^], and with y = m iox right actions and v = n 
for left actions. One can easily prove that: 

• right multiplication MT^ j (p) consists in adding the ith column of M right 
multiplied by p to the jth column of M, the remaining part of M being 
left unchanged, 

• left multiplication Tij{p)M consists in adding the jth line of M left mul- 
tiplied by p to the ith line of M, the remaining part of M being left 
unchanged, 

• TiJ{l)Tj.i{~l)T^.J{^)M (resp. MT^.J{l)Tj.i{-l)TiJ{^)) is the permuta- 
tion matrix replacing the jth line of M by the ith one and replacing the 
jth one of M by the ith one multiplied by —1, all other lines remaining 
unchanged (resp. the permutation matrix replacing the ith cohuim of M 
by the jth one multiplied by —1 and replacing the jth one by the ith one, 
all other columns remaining unchanged). 

Every unimodular matrix V (left) and U (right) may be obtained as a prod- 
uct of such elementary unimodular matrices, possibly with a diagonal matrix 
D{a) = diagjofi, . . . , a.^) with aj G a, ^ 0, i = 1, . . . , z/, at the end since 
T,,,.(p)D(a) = I)(a)r,,,.(^pa,.). 

In addition, every unimodular matrix U is obtained by such a product: 
its decomposition yields VU = I with V finite product of the Ti,j{py& and 
a diagonal matrix. Thus, since the inverse of any Tij{p) is of the same form, 
namely Tij{—p), and since the inverse of a diagonal matrix is diagonal, it results 
that = t/ is a product of elementary matrices of the same form, which 

proves the assertion. 
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B. The Smith-Jacobson Decomposition Algorithm 

The Smith-Jacobson decomposition algorithm of the matrix M consists first 
in permuting columns (resp. lines) to put the element of lowest degree in upper 
left position, denoted by mi^i, or creating this element by euclidean division (in 
A{5)[-^]) of two or more elements of the first line (rcsp. column) by suitable 
right actions (resp. left action). Then right divide all the other elements m\^k 
(resp. left divide the mk,i) of the new first line (resp. first column) by m\^\. 
If one of the rests is non zero, say ri^k (resp. r^^i), subtract the correspond- 
ing column (resp. line) to the first column (resp. line) right multiplied (resp. 
left) by the corresponding quotient q\^k defined by the right euclidean division 
mi^k = TOi,i9i,/c + ''i,fe (resp. qk,i defined by mfe,i = qk^mis+rk^i). Then right 
multiplying all the columns by the corresponding quotients gi,^, k = 2,...,u 
(resp. left multiplying lines by qk^i, k — 2, . . . ,/x), we iterate this process with 
the transformed first line (resp. first column) until it becomes (mi.i, 0, ... ,0) 
(resp. (m,i,i, 0, . . . , 0)"^ where means transposition). We then apply the same 
algorithm to the second line starting from m2.2 and so on. To each transfor- 
mation of lines and columns correspond a left or right elementary unimodular 
matrix and the unimodular matrix V (resp. U) is finally obtained as the product 
of all left (resp. right) elementary unimodular matrices so constructed. 
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